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Magnesium-rich Ejecta in the SNR G284.3—1.8 Around the High-Mass 
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ABSTRACT 

We present results from two Chandra observations of the 16.6 d X-ray/ 7 -ray high- 
mass binary 1FGL J1018.6—5856 located at the center of the supernova remnant 
G284.3—1.8. The binary spectra, separated by 0.25 in binary phase, are fit with an 
absorbed power-law model with T 1.7 — 1.8 for both observations (the flux during 
the second observation is a factor of 1.7 smaller). In the high-resolution ACIS-I image 
we found a hint of extended emission ps 2"-3" southeast of the binary, significant at 
the 3cr level. Binary evolution codes reproduce the system’s observed properties with 
two massive stars with an initial 18-day period, undergoing mass transfer and leaving 
behind a heavy pa 2 M 0 neutron star. The initial mass of the progenitor star in this 
scenario is 27 ±4 M 0 . Chandra and XMM-Newton images of the remnant show it 
has a relatively low X-ray surface brightness. The two brightest regions of extended 
X-ray emission, with luminosities ~ 10 33 erg s _1 for d = 5 kpc, lie in the northern 
and western portions and show significantly different spectra. The northern patch is 
consistent with shocked ISM, with a low temperature and long ionization timescale. 
However, the western patch is dominated by ejecta, and shows significantly enhanced 
Mg content relative to other ejecta products. The abundance ratios inferred resemble 
those from the Large Magellanic Cloud remnant N49B. To our knowledge, this is only 
the second case of such Mg-rich ejecta found in an SNR. Nucleosynthesis models for 
core-collapse SNe predict Mg-rich ejecta from very massive progenitors of > 25 M 0 . 

Subject headings: ISM: supernova remnants — ISM: individual objects (G284.3-1.8) 
— stars: individual (1FGL 1018.6-5856) — X-rays: binaries 
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1. Introduction 


The ejecta from supernova (SN) explosions enrich the interstellar medium of the Galaxy with 
heavy elements over thousands of cubic parsecs. Identification of the ejecta products of SNe is most 
easily done via analysis of the X-ray emission from young to middle-aged supernova remnants 
(SNRs), after the reverse shock has propagated back into the ejecta and heated it to X-ray emitting 
temperatures. One of the biggest open questions in SN research is what types of stellar progenitor 
systems lead to what types of SNe and SNRs. In only a handful of cases have astronomers directly 
observed the pre-SN star in archival images, and this situation is unlikely to change in the near 
future. The fundamental problem is that SNe are plentiful but distant and point-like, while SNRs 
are local and resolvable but limited in number. 


SNR G284.3—1.8 (also kn own as MSH 


dio (Mi lne et al.lll989t) . optical (IRuiz & M ay 


0-53; hereafter G284) has been identified through ra- 


1986), and X-ray (lAbramowski et al.ll2012h observa¬ 


tions, and is perhaps 10 4 years old dRuiz & Mav 1986 ) with a distance of 5.4+ti kpc ( Napoli et al. 


20111) . It is interacting with a molecular cloud in the north. An analysis of the extended X-ray 
emission seen in G284 was first reported by Abramowski et al. (2012), who noted that X-ray emis¬ 
sion was present in the northern portion of the remnant, just behind a filament identifiable in the 
radio and optical images. 

At the center of the remnant lies the bright X-ray point source 1FGL J1018.6—5856 (hereafter 
J1018). Identified as a high-mass binary by Corbet et a l, 2011C J1018 reveals a 16.6 d periodic 
modulation in t he Fermi LAT observat i ons (Ackermann et all2012 ). This is one of the few binaries 
seen in 7 -rays (IKargaltsev et al.ll2014l) . IPavlov et al.1 (1201 lh reported results from the first Chandra 
and XMM-Newton observations of J1018 around binary phases 0 = 0.3 and 0 = 0.6, respectively 
(zero phase corresponds to the maximum of the 1 — 10 GeV flux, MJD 55403.3 ± 0.4), and also 


found the h igh-mass companion in 2MASS data. Suzaku (ITanaka et al.ll2014l) and Swift (lAn et al. 


2013L 2015) observed J1018 covering the entire orbital period, confirming the periodicity seen in 
7 -rays. The most recent period measurement, 16.544 ± 0.008 d, was provided by An et al. (2015), 
who found a large X-ray flux increase at <f = 0 and an additional maximum at f = 0.3 — 0.4. J1018 
exhibits a correlation between X-ray flux and photon index, common to several 7 -ray binaries. 

J1018 is one of two high-mass 7 -ray binaries seen inside of an SNR (along with SS 433 in 
W50), offering a rare opportunity to study the remains of the explosion that would produce such 
a system. The spectral type of the companion star places tighter constraints on the exploded star 
than would otherwise be possible. Ground-based observations revealed a massive star of spectral 
type 06V((f)) (V pa 13.5; Napoli et al. 2011 ) and mass of ~ 30 M 0 , very si milar to that of the 


ray binary LS 5039. J1018 was observed with the H.E.S.S. observatory by Abramowski et al. 


(20150, who report a 9 a detection of the source HESS J1018—589A (spatially coincident with 
J1018). The detected 7 -ray variability, both in GeV and TeV, confirms the association between the 
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H.E.S.S. source and J1018. lAbramowski et ali (12015I) . based on the shape of the TeV phaseogram 
and source spectrum, suggest that J1018 is a low-inclination, low-eccentricity system. 


In this Letter , we report our analysis of the X-ray emission from J1018 and G284. We use 
newly-obtained Chandra data, supplemented with archival XMM-Newton data, to examine the 
extended X-ray emission from the SNR and do spectral and spatial analysis of the binary. 


2. Observations and Data Reduction 


We use two sets of Chandra observations (Table 1), totaling 72 ks, taken with the ACIS-I 
instrument in timed exposure mode and telemetered in “very faint” format. We processed the data 
using the Chandra Interactive Analysis of Observations (CIAcQ) software (ver. 4.6) and Chandra 
Calibration Data Base (CALDB), ver. 4.6.3. Spectra were analyzed using XSPEC (ver. 12.8.2), 
which contains the latest atomic data from AtomDB 3.0. 

We simulate the Chandra PSF using ChaRI@ (Chandra Ray Trace) to study the morphology 
of the emission in the immediate vicinity of J1018. The ChaRT output was supplied to MARX 
v.5.10 software to obtain the final image. We used a blurring (AspectBlur parameter in MARX) 
value of 0719, as recommended by the Chandra Help Desk. 


For our analysis of the G284 emission, discussed below in Section 13.21 we supplement our 
Chandra observations with archival XMM-Newton observations to increase the signal-to-noise ra¬ 
tio. G284 was observed by XMM-Newton on 2013 Jan 9 for a total of 105 ks (ObsID 0694390101; 
PI: De Fuca). EPIC-MOS observations were taken in Full Frame mode with a medium filter, while 
EPIC-pn data were taken in Small Window mode with a thin filter. The pn data captures only emis¬ 
sion from the binary (analyzed by An et al.1120151) . not the extended emission. We processed the 
MOS data with XMM Science Analysis Software (SAS) ver. 13.5. To filter the data for flares, we 
used the mos-filter routine as part of the ESAS software, available within SAS. After time filtering, 
approximately 79 and 91 ks of “good” data remained for the MOS 1 and 2 detectors, respectively. 


Guided by the Chandra images, we selected two regions of extended emission to analyze, 
labeled “north” and “west.” Both are shown in Figure [2l These regions have areas of « 10.9 
and « 16.2 arcmin * 2 , respectively. Background subtraction for these two regions is difficult due 
to the large size of the remnant compared with the Chandra field of view. However, based on a 


'http://cxc.harvard.edu/ciao/index.html 

2 http://cxc.harvard.edu/chart/ 
'http://space.mit.edu/CXC/marx/ 
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comparison of the Chandra and XMM images, we found off-source regions that appear on both the 
ACIS-I and MOS arrays and do not appear to have any associated emission from the SNR. 


3. Results & Discussion 
3.1. 1FGL J1018.6—5856 

3.1.1. Spectral and Imaging Analysis 

We extracted spectra for J1018 from the two Chandra observations. The spectra suffer from 
moderate pile-up (3,294 and 1,498 counts in the r = 3"2 aperture correspond to 0.24 and 0.17 
counts per frame in ObsID 14657 and 16560, respectively). To correct for pile-up, we used the 
Davis (2001) model with the grade migration parameter a = 0.5. Fits with an absorbed power- 
law model (with the XSPEC phabs model) give hydrogen column densities Ah = 9.0 ± 0.9 and 
(9.2 ± 1.1) x 10 21 cm~ 2 ; photon indices T = 1.71 ± 0.12 and 1.77 ± 0.17; absorbed (unabsorbed) 
fluxes To. 5— 8 keV = 1.5 ± 0.1 (2.3 ± 0.2) and 0.9 ± 0.1 (1.4 ± 0.2) x 10 ~ 12 erg s ' 1 cm -2 ; and 
reduced x 2 — 107 and 0.75 for ObsIDs 14657 and 16560, respectively. All uncertainties are la. 
The derived Ah and T do not change significantly between the individual fits. The flux in the 
second observation is a factor of 1.7 smaller than the first one, in accordance with the orbital phase 
dependance of the flux. 

Comparing our fit for the first Chandra observation with the An et al. (2015) fit for the XMM- 
Newton observation taken at a close binary phase (</>cxo = 0.42, </>xmm = 0.33), we find that our 
A h is somewhat larger and the slope is somewhat steeper, but the values are consistent at the 3a 
level. Because of the lack of pile-up, the XMM-Newton measurements are probably more accurate. 

We also investigate the spatial morphology of J1018 in the ACIS-I images and compare the 
surface brightness radial profile with a modeled PSF0 (Figure H). Overall, the observed profile 
matches the simulated PSF. There appears to be a mismatch between the simulated profile and the 
observed profile at r < 1", but this is likely due to an inacurate pile-up modeling in MARX@. In 
addition, the PSF model is not very reliable at small distances due to the presence of the asymmetric 
PSF anomal\@ around r = 078, which is not simulated by MARX. Including the pile-up correction 
in the MARX 5.1 simulation does not change the profile substantially beyond 1", which is expected 


4 We used MARX 5.1 with AspecBlur = 0.19 and pile-up parameter a = 0.5. 
5 http://space. mit.edu/ASCAIARX/indetail/pileup.html 
6 http://cxc. harvard.edu/proposer/POG/html/chap4.html 
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for moderate pile up. 

Analyzing the high-resolution images of J1018, we found an enhancement (24 photons in 
0.5-8 keV within the ellipse shown in Figure 1, for the merged image) located ^ 2"2 away from 
the best-fit centroid of the source at a position angle (PA) of ~ 140° (counted East of North). 
The enhancement, visible in both observations, has ~ 3cr significance estimated from comparison 
with the background (which includes PSF wings) measured in the 1'.'75 < r < 3" annulus for the 
merged image. The expected enhancement due to the known mirror artifact has a very similar 
PA, but is located significantly closer to the source center (r 078) and is unlikely the cause of 
the enhancement. Its appearance is similar to the fast-moving extended feature we recently dis¬ 
covered in the vicinity of B1259-63 7 -ray binary (Pavlov et al. 2015). Deeper followup Chandra 
observations are needed to study this extended emission and its tentative connection to J1018. 


3.1.2. Modeling of Binary Evolution 


Although the optical counterpart of J1018 has been found and the binary nature of this system 
is firmly established, it remains unknown whether the compact object is a neutron star (NS) or a 
black hole (BH). We use the known binary properties (orbital period and mass of the donor star) 
to place constraints on the compact object. While there is still some uncertainty about the exact 
range of masses of progenitors that end their lives as NSs and those that become BHs, most models 
predict that, for solar or lower metallicities, the transition between NSs and BHs occurs for stars 
with initial masses somewhere in the range of ~ 18 — 25 M 0 (e.g.. [Fryer 1999; Heger et a,l. 20031) . 
For metallicities much larger than solar, stellar evolution models predict that massive stars develop 
substantial winds, which cause enough mass loss that the end product is a NS rather than a BH 
(Heger et al. 20031) . For our simulation we assume solar metallicity as the system is in the Galactic 
disk. We first use the single star evolution (SSE) code (Hurley et al. 2000) to study the evolution of 
stars with solar metallicity. We find that stars with initial mass of 20 M e and higher would end up 
as BHs, in agreement with the above. If J1018 started and evolved as a detached system (no mass 
transfer ever happens), the progenitor star must have burned its hydrogen faster than the 3OM 0 star 
and hence it should have been even more massive. Therefore, in this scenario the compact object 
should be a BH. However, the lack of mass transfer would require a much larger than the currently 
observed (16.6 d) initial orbital period, which could only increase (since the binary is assumed 
wide and there is no interaction to cause spiraling in) after the progenitor explodes as supernova. 
Thus we conclude that the system must have started as a tight binary and undergone mass transfer. 


7 http://cxc. harvard.edu/ciao/caveats/psf_artifact.html 
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However, evolution of massive stars in tight binaries is more complicated because mass trans¬ 
fer can proceed from one star to the other, rendering the SSE code inapplicable. Therefore, we 
used binary stellar evolution (BSE) models (Hurley et al. 2002). We ran a grid of 80,000 simula¬ 
tions with the following parameters: M\ = 10 — 35 M 0 (current donor star), M 2 = 12 — 50 M 0 
(SN/compact object progenitor with initial mass always greater than Ml), binary period P = 
5 — 50 d, eccentricity e = 0 — 0.9. We find 115 systems that end up with P = 16.6 ± 1.5 d, 
and Mi = 31 ± 3 M 0 right after the SN explosion of star 2. A large fraction of these simulated 
binaries are disrupted immediately after the SN, but we are not interested in such systems. Our 
findings show that in order to reproduce the observed binary parameters we need the following 
initial conditions: Mj ni = 13.4 ± 2.5 M 0 , M™ 1 = 26.7 ± 3.7 M 0 , P = 18 ± 11 d. The simulations 
resulted in a heavy (2.2 ± 0.4) NS as the compact object in 95% of the cases, while the 06V donor 
star gained mass (from 13.4 M 0 to 31 M 0 ) via mass transfer. Note that the sys tems that match 
the above criteria have high eccentricity (e = 0.57 ± 0.24), while lAbramowski et ali (|2015 ) claim 
(based on the TeV observations) that the eccentricity of J1018 should be closer to 0. 


3.2. SNR G284.3—1.8 


The northern region analyzed in Abramowski et al. (2012) roughly corresponds to our “north” 
region. We extend their analysis, supplementing it with additional Chandra and deeper XMM- 
Newton data (their analysis used a 20 ks XMM observation), while expanding the analysis to our 
“west” region as well. While we do not show it here, we have also found that the northern filament, 
identifiable in radio, Ha, and X-rays, is present in the WISE 22 pm IR data. IR emission from the 
forward shock in SNRs is common, even in more evolved SNRs ( Sankrit et all2014 ). and generally 
results from dust heated in the post-shock environment. 

We extracted spectra from the north region using only the MOS 2 and Chandra observations, 
since this region falls on one of the dead MOS 1 chips. For the west region, data from MOS 1, 
MOS 2, and Chandra are used. After background subtraction, slight residuals remained in the 
MOS 1 & 2 spectra from instrumental lines of Al and Si Ka at 1.49 and 1.75 keV, respectively. 
We excised small regions around both of these lines (1.45-1.53 keV and 1.71-1.79 keV) from the 
MOS 1 & 2 spectra, but kept those excised regions in the Chandra spectra. For each region, we 
perform joint fits to all available data, with all parameters tied together. The column density, is 
taken from the fits to the X-ray spectrum of the point source, and is frozen in the fits to 8xl0 21 


cm , the same value as lAbramowski et all |2012) used. We binned the data to a minimum of 50 


counts per spectral bin for all SNR datasets. 


We show the spectra from the north and west regions in Figure [3} We fit these spectra with 
an absorbed plane shock model with variable abundances (phabs*vpshock), assuming cosmic 
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abundances from Wilms et al. ( 20021) . The temperature and ionization timescale (r = j'J n e dt ) 
were allowed to float freely, as were abundances of Ne, Mg, and Fe. Si was free as well, but was 
tied together with S due to the poor statistics of the S line. We kept O frozen at the solar value for 
two reasons: with the fairly high absorption towards the remnant, O lines are only barely visible, 
if present at all, and we wanted the abundances of the other elements to be normalized to some 
fiducial value (the abundance of O). 


For the north region, we find normal cosmic abundances for Ne and Mg, with subsolar abun¬ 
dances for Si (TS) and Fe. Subsolar abundances are often seen in the X-ray spectra of middle-aged 
SNRs such as the Cygnus Loop ( Tsunemi et al.ll2007 ). T he relatively long ionizatio n timescale and 
lower temperature leads us to concur with the analysis of Abramowski et al. (2012) that this emis¬ 
sion is dominated by shocked ISM. The ionization timescale is fit as 4.6 x 10 12 cm -3 s, but is 
unbound on the upper end (lower-bound of 1.9 xlO 12 ). An ionization timescale this long makes 
the spectrum almost indistinguishable from a plasma in collisional ionization equilibrium, again, 
not unexpected for the shocked ISM from an older remnant. 


The spectra of the west region show an entirely different emission profile. They are relatively 
well-fit (reduced x 2 of 1.4 for 184 d.o.f.) by the same absorbed plane-shock model with variable 
abundances with the parameters shown in Table 0 Immediately obvious is the strong Mg line 
at 1.35 keV. Ne and Fe are both consistent within errors with solar abundances, while Si (and S) 
may be slightly overabundant. However, the overabundance of Mg is about a factor of 4.5 higher 
with respect to O. This overabundance, combined with the higher temperature (~ 1 keV) and 
much lower ionization timescale (~ 1 x 10 11 cm ~ 3 s) implies that the emission in this region is 
dominated by Mg-rich ejecta, heated by the reverse shock. This conclusion is not sensitive to the 
choice of model; a regular NEI model with variable abundances (model vnei in XSpec) returns 
similar values. 


The spectral fits for this region are similar to those reported in Park et al.1 (2 0031) for the LMC 
remnant N49B. Their analysis is similar to ours, albeit with LMC abundances, different from 
those in the Galaxy. They find an enrichment of Mg and Si (though their Si abundance is higher 
than ours), while Ne and Fe are not enhanced. They measure a Mg/O ratio of 4.2 from a small 
portion in the interior of N49B, within errors of our measurement of 4.5 in G284. As pointed 
out by Park et al. ( 20031) . the high Mg/O ratio is surprising, as nucleosynthetic models predict that 
significant quantities of Mg produced in core-collapse SNe should be accompanied by significant 
quantities of O ( Thielemann et al.lll996 ). While examples of significant O without Mg are plentiful, 
it is unexpected to find significant Mg without O. 


As Park et al. (2 003b show, deriving a total mass of Mg ejected by the SN from X-ray spectral 
analysis is very difficult, with estimates varying by more than an order of magnitude depending on 
the assumptions made. Such an estimate is beyond the scope of this work. Nonetheless, the high 



















- 8 - 


Mg abundances in the ejecta suggest a substantial amount was synthesized in the explosion. In the 
nucleosynthesis models of Thielemann et al. ( 19961) , significant amounts of Mg are only produced 
in massive progenitors of > 25 M 0 . 


4. Conclusions 

While the northern region of G284 is dominated by shocked ISM, we have shown that the 
ejecta-dominated western region is significantly overabundant in Mg. Along with N49B in the 
LMC, this makes only the second known remnant with Mg-rich ejecta. Models of nucleosynthesis 
in CCSN progenitors produce significant amounts of Mg relative to other products, like Si and Fe, 
only in the explosions of massive stars (> 25 M 0 ). 

We find no statistically significant extended emission in the radial profile for r < 7" from the 
binary J1018 at the center of G284, though we do see a slight excess of photons a few arcseconds 
southeast of the source. The spectral parameters of J1018 are about the same in the two observa¬ 
tions shifted by 0.25 in binary phase, but the flux difference is significant, « 30-40%. Our binary 
evolution calculations support the idea that the progenitor for this SN had an initial mass of > 25 
M q . Mass transfer prior to the explosion from the progenitor to the current donor star increased 
that star’s mass to ~ 30 M 0 . Our simulations imply that the compact object in the binary is a 
heavy, ~ 2M 0 , NS. 

With its large size and low surface brightness, G284 would benefit from a deep study with an 
X-ray telescope with both large effective area and field of view. A long observation with the pn 
instrument on XMM-Newton in Full Frame mode or the Soft X-ray Imager (SXI) on the upcoming 
Astro-H satellite would be ideal for the purpose of further study on the fainter sections of the ejecta. 


We thank the anonymous referee for useful comments. We acknowledge useful discussions 
with Steve Snowden on the XMM-Newton Extended Source Analysis Software (ESAS), a part 
of the SAS software distributed by the XMM-Newton Guest Observer Facility. We acknowledge 
Chandra Grant G03-14047. 
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Table 1. List of Chandra observations 


Obsld 

Date 

</> a 

Exp b 

14657 

2013-12-24 

0.424 ± 0.044 

44.40(0.031) 

16560 

2013-12-29 

0.671 ± 0.044 

27.68(0.019) 


a Orbital phase of the middle of the observational in¬ 
terval. Zero phase corresponds to the 7 -ray peak, MJD 
55403.3_± 0.04; the binary period is 16.544 ± 0.008 d 
( An et al. 2015 ) 


b Exposure in units of ks, and the corresponding phase 
interval (in parentheses). 
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Table 2. X-ray Spectral Fits for G284 


Region 

kT (keV) 

Ne 

Mg 

Si (S) 

Fe 

n (cm 3 s) 

Norm. 

Flux 

xl (d.o.f.) 

North 

0.67 3 ; 6 1 

1 IQ 1 - 70 
a - a *0.86 

1-06 bf 9 

0 IQ 0 - 32 
u * iy 0.09 

0.24°; 3 1 

4 . 61.9 xlO 12 

9.6\ 1 j 7 xlO " 4 

3.1 

1.96(112) 

West 

0.924;°? 

1 30 1 - 72 
a * ju 0.95 

A cq6.39 
^• JJ 3A9 

i cn2.16 
a * ju 1.08 

0-97 2 ;i 

l.Ogi xlO 11 

1.6g;f xlO " 4 

2.5 

1.40(184) 


Note. — Fits to the two regions of G284 shown in Figured Fits were performed using XSPEC with a phabs * vpshock 
model. Abundances are relative to solar. Subscripts and supercripts are lower and upper bounds, respectively, at the 90% 
confidence level. Flux in units of 10 -13 erg s - ' 1 cm ~ 2 over range 0.5 — 2.5 keV. 
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Fig. 1.— Left: ACIS-I merged image of J1018 (pixel size 0"25). Small black circle shows the best- 
fit centroid position. The black sectors show the location of the known Chandra mirror artifact. A 
slight excess is seen within the white ellipse, see text for details. North is up, East is to the left. 
Right: Observed surface brightness radial profile of J1018 (black) and modeled PSF (red) with 
account for pile-up, see Section 3.1.1 for details. 



Fig. 2.— Left: Chandra image of G284, with 0.5-1.2 keV emission in red, 1.2-2.0 keV in green, 
and 2.0-7.0 keV in blue. Right: XMM-Newton EPIC-MOS 2 image of the remnant, with identical 
color bands. Overlaid are the extraction regions, with the “north” region in cyan, the “west” in 
white, and the backgrounds in green. 
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Fig. 3.— Left: X-ray spectra from the “north” region of G284, with Chandra ACIS-I data shown 
in black and XMM-Newton MOS 2 data shown in red. A fit from a single absorbed vpshock model 
is overlaid. Fits were performed jointly to both data sets simultaneously using XSPEC. Right: 
X-ray spectra from the “west” region of G284, with Chandra ACIS-I data shown in black, and 
XMM-Newton MOS 1 & 2 shown in red and green, respectively. Overlaid is an absorbed vpshock 
model fit. 
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